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The kinetic approach is applied to develop the Drude-Sommerfeld model for studying of the optical 
and electrical transport properties of spheroidal metallic nanoparticles, when the free electron path 
is much greater than the particle size. For the nanoparticles of an oblate or a prolate spheroidal 
shape there has been found the dependence of the dielectric function and the electric conductivity 
on a number of factors, including the frequency, the particle radius, the spheroidal aspect ratio and 
the orientation of the electric field with respect to the particle axes. The oscillations of the real 
and imaginary parts of the dielectric permeability have been found with increasing of particle size 
at some fixed frequencies or with frequency increasing at some fixed radius of a nanoparticle. The 
results obtained in kinetic approach are compared with the known from the classical model. 
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I. INTRODUCTION 

Understanding how light interacts with matter at the 
nanometre scale is a fundamental problem in nanopho- 
tonics and optoelectronics. The variety of applications 
are based on the surface-bound optical excitations in 
metallic nanoparticles (MNs). They include plasmonic 
nanomaterials, optical nanosensors, biomarkers, inte- 
grated circuits and sub- wavelength waveguides, molec- 
ular imaging, metamaterials with a negative refractive 
index, optical antennas and so on (see, e.g., Refs.[l|-[2|). 

The optical and transport properties of the MNs 
have been investigated for a long time and the results 
have been presented quite completely in numbers of 



reviews*— and monographs 



12-16 



especially for MNs of 



a spherical shape. 

Dielectric function is the most important factor for 
the design and optomization of plasmon nanometer-sized 
structures. However, this function for MNs differs from 
an ideal bulk metal. The difference depends on many fac- 
tors, such as the size and shape of the MNs and the sur- 
rounding mediaJ^^— That is why the dielectric function 
of MNs has been under comprehensive study for many 
years.— ^— The size effect becomes more significant when 
the size d of MN is comparable with the electron mean 
free path I. The ratio between d and I is a very impor- 
tant physical characteristic, which eventually defines the 
mechanism of an electron relaxation rate inside the MN. 
As a rule, the cases d ^ I and/or d <C I are studied. The 
former case refers to the so called diffusive electron dy- 
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namics and the latter one to the ballistic electron dynam- 
ics. The diffusive case has been studied in details, since 
it enables to use the Mie theory for a uniform media3^ or 
the Maxwell- Garnett theorySl for a composite media in 
the calculations of the optical properties of the MN. 

In particular, when d ^ I, the following expression for 
the dielectric permeability which results from the Drude- 
Sommerfeld theoryi^iS^ has been often used: 



e(w)==e'(w)+ie"(w) = l- 
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Here v is the collision frequency inside the particle bulk, 
LUpi = ^J Anne^ /m is the frequency of plasma electrons 
oscillations in the metal, e and m are the electron charge 
and mass, correspondingly, and n is the electron concen- 
tration. The imaginary part of the dielectric permeabil- 
ity e"(w) is connected with the high-frequency (optical) 
conductivity by the well known relation: 



a{.) = -e"{u:) 
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(2) 



For sufficiently low frequencies (w -^ 0), one gets the 
expression ctq — ne^ /mv describing the statical conduc- 
tivity. 

In the case when the sizes of the MN are less than the 
mean electronic free path in the particle, the mechanism 
of electrons scattering is changed, and the surface of the 
particle starts to play a dominate role. This effect will be 
called further as the surface or boundary effect. Strictly 
speaking, in this case neither Mie nor Drude-Sommerfeld 
theory can be applied and an another approach ought to 
be elaborated to present the optical properties of MNs. 
As it was shown^i21 a clear size dependence of relaxation 
dynamics is observed in experiments. This is a strong 



indication of an efficient electron-surface phonon interac- 
tion in this regime. But the experimental results^i^Sr^ 
are often analyzed within the frame of the named theo- 
ries, or to simplify the problem, the Eqs. (HJ, 1^ have 
been used with formal replacementsi^i^"— 
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where vp is the electron velocity at the Fermi surface, R 
refers the particle radius, V and S refer the volume and 
the surface area of the spherical particle, and ^ is a coef- 
ficient obtained by fitting the calculations to the exper- 
imental data. In other words, the collision frequency in 
the particle volume is replaced by some effective collision 
frequency. But such a replacement can be applied for the 
MNs of a spherical shape only. If the geometry of MN 
differs from the spherical one, and the condition d <^ I 
holds for at least one of the particle directions, then the 
optical conductivity becomes the tensor quantity^! and 
the formal replacement similar to the one presented by 
Eq. ^ can not be used any more. It is necessary to 
look more closely at the effect of the particle boundaries 
on the optical properties. For a theoretical study it is 
convenient to choose the particle of a spheroidal shape. 
That is because the results obtained for the particles of 
such a shape can be easily extended to the particles of 
other shapes by means of the formal transformation of 
the spheroidal axes. 

Thus, there arises a problem in the case d <^ I, how to 
calculate the values of e or a, when the mentioned above 
classical models can not be used. 

The present work is devoted to the elaboration of the 
new approach, allowing to calculate the real and imag- 
inary parts of the dielectric permeability for MNs of a 
spheroidal shape in the case when the inequality d <C I 
is hold. 

Since both the real and the imaginary parts of e(w) 
govern the numerous properties of the MNs, the neces- 
sity of the detailed study of the boundary effect on the 
dielectric permeability has become evident. 

The rest of the paper is organized as follows. The 
kinetic approach to the problem is presented in Sec. II. 
Section III contains the study of the conductivity in the 
spheroidal MNs. In Sec. IV, we consider the limit cases 
of the problem and the size effects. Sec. V is devoted 
to the discussion of the obtained results, and Sec. VI 
contains the conclusions. 



II. KINETIC EQUATION METHOD 

To account for the effect of MN boundaries on the op- 
tical and electrical transport properties of the MN, we 
will apply the kinetic equations approach. The advan- 
tages of this approach is that the obtained results can be 
applied to strongly anisotropic spheroidal (needle-like or 
disk-shaped) MNs, but in the case of MNs of spherical 
shape it transforms to the well known resultsi^— , like 



the ones given by Eqs. ([I]), ([2]). Thus, it permits one to 
study the effect of the particle shape on the measured 
physical values. Secondly, the kinetic method enables to 
investigate the MNs with sizes as more or less then the 
electron mean- free path I. But, there exists the lower 
limit of the applicability of this method in the small ra- 
dius limit, when the particle size is comparable with the 
de Broglie wavelength of the electron, and the quantiza- 
tion of the electron spectrum start to play an essential 
rolci^"— Practically, it is around a radius of more or less 
2 nm. Beyond the electronic Fermi distribution function, 
the influence of quantum effects on the optical conduc- 
tivity can manifest itself in the quantization of electron 
pulses and angular momentums. It is important at low 
temperatures, when the distance between successive en- 
ergy levels is much higher than ksT. 

Let us consider the single MN which is irradiated by 
electromagnetic wave, whose electric field is given as 



E = Eoexp[i(kr - tot)]. 



(4) 



Here Eq is the amplitude of an electric field, ut is its 
frequency, k is the wave vector, and r and t describe the 
spatial coordinates and time. 

We will assume that the electromagnetic wave length 
is far above the particle size. If one chooses the coordi- 
nate origin in the center of the particle, then the above 
mentioned assumption is written as 



kr < 1. 



(5) 



The inequality ([S]) implies that the E-field of the elec- 
tromagnetic wave can be considered as spatially uniform 
on scales of the order of a particle size. This means that 
the field represented by Eq. (JH) induces inside of the MN 
the electric field varying in time, but uniform in space. 
The amplitude of such a field is connected with Eq by 
the relation"^^ 



ELi^) = 



E'o{0,i^) 



1 + L,[eiuj) - ly 
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where Lj are depolarization factors in the j-th direction 
(in the principal axes of an ellipsoid). The explicit ex- 
pressions of Lj for a single MN of an ellipsoidal shape 
can be found elsewhere (see, e.g., Refs. [40.41]). 

The field Ein has an effect on the equilibrium electron 
velocities distribution and, thus, determines the appear- 
ance of a nonequilibrium addition /i(r, v, t) to the Fermi 
distribution function fo{e). Here e = mv'^/2 is the ki- 
netic energy of electron, and u = |v| refers the electron 
velocity. As is well known,— the equilibrium function 
/0(e) does not give any input to the current. With the 
account for both the time dependence of Eq. ^ and 
the inequality ([5]), the distribution function of electrons, 
which generates the field Ei„, can be written as 



/(r,v,t) = /o(e) + /i(r,v,i) = /o(£) + /i(r,v)e'-*. (7) 



The function /i(r, v) can be found as a solution of the 
Unearized Boltzmann's equation 
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In Eq. ([5]) we have assumed that the colhsion integral 
{dfi/dt)coi = —Ji/t is evaluated in the relaxation time 
approximation (r = 1/^^)- What is more, the function 
/i(r, v) ought to satisfy the boundary conditions as well. 
These conditions may be chosen from the character of 
electrons reflection from the inner walls of the MN. We 
will take, as is usually done, the assumption of a diffusive 
electron scattering by the boundary of MN. Placing the 
origin of the coordinates at the center of the particle, we 
can present this boundary conditions in the form 



/(r,v)|s=0. 



<0, 



(9) 



where v„ is the velocity of the component normal to the 
particle surface. 

Alongside with diffusive, the mirror boundary condi- 
tions at the nanoparticle surface were examined in liter- 
ature for electron scattering (see, e.g., Refs. [42, |4J]). In 
this case, each electron is reflected from the surface at 
the same angle at which it falls to the surface. In diffuse 
reflection, the electron is reflected from the surface at any 
angle. In order that the mirror mechanism be dominant, 
the surface must be perfectly smooth in the atomic scale, 
since the degree of reflectivity of the boundary essentially 
depends on its smoothness. Practically, for a nonplanar 
border such smoothness is extremely difficult to achieve. 
As it was shown**^, the mirror boundary conditions give 
a small corrections to the results obtained with the ac- 
count of only the diffusive electron reflections. There- 
fore, we chose more realistic boundary conditions given 
by Eq. ^. 

It is comparatively easy to solve Eq. ^ and to satisfy 
the boundary conditions of Eq. (|9l), if one passes to the 
transformed coordinate system, where an ellipsoid with 
semiaxes i?i, i?2, ^3 transforms into a sphere of radius R 
with the same total volume: 

^' ' R^{R,R2R,f\ (10) 






Similar transformation should be made for the electron 
velocities as well: Vj — v'.Rj/R. Then solving Eq. ([5]), 
one finds 

/i(.r,v,tj — -e— — vJiin ; , 

(11) 
where the characteristic of Eq. ([5]) can be presented as 



Generally speaking, the presence of the surface changes 
both the current and the field distributions. It is reason- 
able to point out here that, though the electric field still 
remains homogeneous inside of the MN (in accordance 
with the Eq. ([5])), the distribution function /i is the co- 
ordinate dependent in any case due to the necessity to 
meet the boundary condition given by Eq. ^. 

Performing the Fourier transformation of Eq. dill) , one 
can calculate the density of a high-frequency current in- 
duced by the electromagnetic wave of Eq. @ inside the 
MN, via the expression 






v/i(r,v,a;)d v. 



(13) 



Let us introduce the tensor of the complex conductivity 
cr^o(r,aj) using the relationship 



13=1 



(14) 



Then in accordance with the both Eq. (|TT]) and Eq. (|13p , 
the components of this tensor can be presented in the 
form 



aS,(r,-) = 2e(^) 
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(15) 



Before the detailed study of the role of the particle sur- 
face and the size effects, which will be displayed below, 
it is worth to note here the following. The surface effect 
on the conducting phenomenon is described in Eq. (jlSp 
by means of the characteristic i'(r', v'). It accounts for 
the restrictions imposed on the electron movement by a 
nanoparticle surfaces. As one can see from Eq. ([T^ . the 
value of t' is of the order of t' ~ R/vp, where vp is the 
Fermi velocity. This implies that the value reciprocal to 
t' will correspond to the vibration frequency between the 
particle walls. Hence, the inequality vt' 3> 1 indicates 
that the electron collision frequency inside the volume of 
MN would significantly exceed the one for an electron 
collision with the surface of MN. If pointed inequality is 
satisfied, it is possible to direct t' — >■ oo and the expo- 
nent in Eq. ([15]) can be neglected. Then, we obtain a 
standard expression for the dielectric permeability, such 
as given in Eq. ^. To ensure in that, it is necessary to 
pass in Eq. (J15p to the integration over v in the spherical 
coordinate system 



V ^ 



(R2-r'2)v'2 + (r'v')2 



(12) 



The characteristic curve of Eq. (|T2|) depends only on the 
absolute value of R and does not depend on the direc- 
tion of R. The radius vector R determines the starting 
position of an electron at the moment t' = 0. 
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with the use of the following formulas 



^5{v'-vl)dv^^-f, 



(16) 



T V 27rfi 



(17) 



and take into account that the energy derivative of /o in 
zero approximation in the smah ratio of ksT/ep, can be 
replaced by 



^^0 x/ \ 



(18) 



as well as the fact that only the diagonal terms with 
Va — vp — V are retained after integration over all angles. 
At the end of this section, we would like to pay atten- 
tion on the next two important circumstances, i) Though 
the inner field Ei„ is spatially uniform, the distribution 
function /i(r,v,t), and, consequently, the density of the 
current j(r,aj) are depended on the coordinates. This 
dependence is imposed by the boundary conditions of 
Eq. ([8]). ii) The density of the current and the compo- 
nents of the conductivity tensor have the physical sense 
only if they are averaged over the particle volume. For 
instance, it is easy to ensure that the energy absorbed by 
a single MN is determined either by an averaged density 
of the current (j) or by an averaged tensor of the com- 
plex conductivity {a^a{r,uj)) (with the account of the 
Eq. dUD). 



III. CONDUCTIVITY OF SPHEROIDAL 
METALLIC NANOPARTICLE 

Let us average over coordinates the components of 
the conductivity complex tensor represented by Eq. (ITSI) . 
The necessity of such an averaging arises from the fact 
that the power absorbed by a single MN is caused by the 
conductivity averaged over the volume of MN. Then, in 
view of Eq. P^ . one gets the expression 



X J d^VVaVp Slv"^ — Vp) 1 
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where the angle brackets denotes the averaging. Firstly, 
we can fulfil the integration in Eq. (jTH]) over all electron 
coordinates. In accordance with Ref. [39| , one can find 
that 
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where the following notations have been used 
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{V - luj) = 91 - iq2- 



(21) 



(22) 



Further, we will take into account only the diagonal com- 
ponents of the conductivity tensor. Based on Eq. (I^Ul) . 



we can rewrite Eq. p9)) as 

Ka(l-,W)) - 
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(23) 



Let us restrict ourselves only with the single MN of 
spheroidal shape. To calculate the residual integral, we 
pass to the spherical coordinate system with a z-axis di- 
rected along the rotation axis of the spheroid (as we have 
done it in the previous section), and take into account 
that the components of an electron velocity parallel (un) 
and perpendicular (uj_) to this axis are defined as 



wii = Vz = V cost 



v± =^ Jvl + v^^v sin e, (24) 



correspondingly, where 
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After integration in Eq. (I23p over the azimuthal angle Lp 
and over all electron velocities, in view of Eqs. p^ . (|19p . 
we obtain for the ||- and _L-components of the conductiv- 
ity tensor the expressions 

C^ii = K;.(l-,W)) =7''^ 
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(25) 



and 
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The subscript v = vp means that the electron velocity 
in the final expressions should be taken on the Fermi 
surface. The ^-function in Eqs. (pS)) . ([25)) varies now 
with the angle 9, because q (see Eq. ([H])) for a spheroidal 
particle becomes dependent on the angle 9, and can be 
determined as 
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where Rn and R^ are the semiaxes of the spheroid. Such 
form for q follows from the form of " deformed" electron 
velocity, which enters into Eq. ([^^ , and for a spheroid is 



v' — vR\ 
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v'{e), (28) 



where the velocity components f ii and v_\_, presented by 
Eq. (|M| . were used. The velocity v' does not depend on 



the particle radius R, but only on the spheroid aspect The calculations were carried out using such parameters 



ratio. In the case of a spherical particle i?|| — R± = R, 
and v' — V. 

There is a well known common relation'*^ between the 
tensor components of the complex dielectric permeability 
and the components of a complex conductivity tensor 

47r 
{ecfi{r,uj)) =Scfi + i—{cr''^p{r,uj)). (29) 

If one separates the real and the imaginary parts in the 
expressions for both the complex tensors of the dielectric 
permeability and the conductivity, i.e., presents 



for Au: v = 3.39 x 10" at 0°C [39,], n = 5.9 x lO^^^-^-s 
VF = 1.39 X lO^cm/s ^. 
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then in a correspondence with Eq. (p9|). one obtains the 
next two relations 



and 
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Finally, using Eqs. (|25]), ^ and Eq. ([31]), one gets for 
a spheroidal particle 
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(35) 
The upper (lower) symbols in the parentheses of the left 
hand sides of Eqs. ([M]) . (1551) correspond to the upper 
(lower) symbols in the parentheses of the right hand sides 
of these equations. Formulas ([M]). ([55]) are the funda- 
mental equations for calculations of the dielectric func- 
tion and, therefore, for studying the optical properties of 
MNs. 

For illustration, we present in Fig. 1 the frequency 
dependence of an imaginary part of the dielectric per- 
meability ratio components for a spheroidal Au particle, 
which is obtained by numerical evaluating the integrals 
in Eqs. ([55)) . ([M)) . It is worth to note that the radius i?|| 
is directed along the revolution axis of the spheroid, and 
R± - transverse to it. These spheroidal radiuses can be 
expressed through the radius of a sphere (of an equivava- 
lent volume) as 

R^{R.Rl)^/\ 




FIG. 1. The dependence of the ratio ^'\_h'a for prolate (_Rx/^|| = 
0.1, upper curve) and oblate {R_i_/Rn = 10, lower curve) Au 
particle (with R = 50A) versus frequencies ratio lo/us, where 
Us = vp/{2R). 

As one can see, the ratio of fj^/eif oscillates with in- 
creasing frequency both for the prolate and the oblate Au 
nanoparticles. These oscillations have a damping char- 
acter both for prolate and oblate particles, but differ in 
the period of oscillations. For a given prolate nanopar- 
ticle the oscillations occur around the constant value 
e" /ejf ~ 4/3 and for an oblate ones - in the vicinity of the 

another constant f"^/f'u ~ 1/2. The period of oscillations 
depend only slightly on the particle volume, however, the 
oscillation amplitude is more pronounced for particles of 
smaller radii. 

It is also worth noting that the intensity of the surface 
mode is determined by the magnitude of the imaginary 
component of the material dielectric constant. Materials 
with a small e" have a large, narrow absorption peak, 
whereas materials with a large e" have a small broader 
absorption peak<^ 

The ratio of the real parts of e'^/eji does not oscillate 
with frequency, therewith at the frequency of w ~ ujpi 
exhibits the singularity. Thus, we demonstrate below the 
plots for e^(a;) and e'Juj) separately. 

Figure 2 shows the real part of the dielectric perme- 
ability components for Au nanoparticle as a function of 
normalized frequency, obtained by using Eqs. ([5^ . ([55)) 
in numerical calculations. The magnitude of e' reaches 
minimum value at a; — >■ 0, and e' — ^ at w — >■ Wp/. Among 
the two components of e', the frequency dependence is 
more pronounced for the longitudinal one: en attains the 
smallest negative value at cj — > in the case of the pro- 
late Au nanoparticle, and the maximal negative value - 
for the oblate one. The absolute magnitude of the both 
components of e' is essentially enhanced as the radius of 
the particle is increased (especially, at w — )• 0). 
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FIG. 2. The dependence of £^'^11 f°^ both the prolate (Rx/^|| = 
0.1, sohd lines) and the oblate {R±/Ru = 10, dashed lines) Au par- 
ticle (with R = 50A) versus frequencies ratio uj/us- Thick curves 
are for ||-components, and thin curves - for ±-coinponents of e' . 

Below, we will consider several approaches, enable us 
to derive the explicit analytical expressions for a' and a" 
froniEqs. dMl, (|33- 



IV. LIMIT CASES. SIZE EFFECTS 



which enters into Eqs. (p4)) . ((35|) . Here, the velocity v'{0) 
is given by Eq. (gS]). The formula dSHl) at z^ -> agrees 
with an earlier estimation given in Ref. p5l |. 

Substituting Eq. §^ into Eq. dM]), and Eq. ^ into 
Eq. psp , correspondingly, and using the values of the in- 



tegrals /|| and I± given in the Appendix A, one obtains 
for the real and imaginary parts of the conductivity ten- 
sor components the expressions 
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provided that the conditions v ^ Vg, uj ^ Vs and v — vp 
are satisfied. To meet the requirement of j^ ^ I's, e.g., 
for Au particle at 0°C, it is necessary that its radius 
should be i? ^ 400A. The first term in the parentheses 
of Eqs. (|^T|) . (|^^ results from the integration over angle 



-it/2 



fsm cos^ 
isin^ 61 



de 



(43) 



A. Frequency approach 

The straightforward evaluations of the dielectric per- 
meability of a single MN or it conductance can be made, 
when ^-function entering in Eq. (1^^ or in Eqs. (IMl) . ([55)1 
takes the simplest form. Let us introduce the value 



V 

2R' 



(36) 



which will characterize the frequency of electron collisions 
with the spherical MN surfaces, and R is the sphere ra- 
dius. This allow us to rewrite Eq. (22) as 



Vs 



(37) 



i) First, we consider, for example, the case |g| ^ 1. 
In the frequency scale, this implies that both inequalities 
V ':^ Vs and uj '^ Vs must be executed. Then, Eq. (PT|) 
reduces to the form 
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Accounting for Eqs. (|3S|) and ([22|) - we can calculate ap- 
proximately the real and imaginary parts of the ratio 
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The second term in the parentheses of Eqs. (|4T|) . (j42|) 
contains the integrals / which accounts for the particle 
nonsphericity and doesn't depend on the frequency. The 
simplest result for a can be obtained for a spherical MN, 
when the integrals / for different polarizations coincide 
with each other, and are equal to 1/3 (see Eqs. (jASp in 
Appendix A). 

As one can see, the first term in the parentheses of 
Eqs. (HI]), (J32I) describes the Drude-Sommerfeld results 
for a spherical particle, and the second one gives the 
first correction of the kinetic theory to the volume elec- 
tron scattering, allowing to account an electron scatter- 
ing from the surfaces of the nanoparticle as well. The 
calculations of the frequency dependencies of the ra- 
tio e'i/eJi with employing of Eqs. ([22]) and (jH]) (and 
the conditions pointed therein) give the results qualita- 
tively similar to our previous numerical calculations (us- 
ing Eqs. (P51) . ([M]) ). depicted in Fig. 1, though without 
any oscillations. But quantitatively, the results obtained 
from Eq. p4l) are of one or even two orders of magnitude 
lower (depending on the radius of MN) than those fol- 
lowing else from the Drude-Sommerfeld formula or from 
Eqs. dH]). 

ii) For the opposite limit case, when \q\ <C 1 (or v <C Vs, 
and oj <C i^s), we can take advantage of the expansion 
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'L 
3! 
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4! 



5! 



6! 



and then from Eq. (1211) . we immediately obtain 
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In this case, we can find for real and imaginary parts of 
the mentioned above ratio the relations 



^ 



^(9) 



R 
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, (45) 



(46) 



This permits us to get for a components the following 
expressions: 



^9ne2 fR 2 y 



(1) 



6\) m vi \ \-i-) 



(47) 



(48) 



provided that i^ <C Vs, uj <C i^s, and v — vf- In other 
words, the first condition means that the radius of MN 
should obey the condition R ^ vp/i'^v)- The first term 
in the parentheses of Eqs. (|47|) is the result of the inte- 
gration over angle 9 with v'{6) in the denominator (see 
Appendix A, Eqs. ()A6p . (jA7|) ). and the second one can 
be given by 
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\{R^IR\\Y 



arctan 
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R\-Rl 
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(49) 



The upper (lower) signs and upper (lower) symbols in the 
parentheses of the right hand side of Eq. (|49p correspond 
to those in the left hand side of this equation. In the case 
of prolate spheroidal particles (i?|| > R\_)^ one should 
made in Eqs. (H^ only the following replacement 



arctan 



i?' -Hi 



R\ 



2i 



R\\ 



-ri|| — li^ 



R\\ 



Rj-Rl 



, (50) 



'(!) 



and in the case of a spherical symmetric MN, _R|| = Ri 
R, and J|| ^ J± = 1/3. 

As one can see from Eq. (|47p . the real part of cr 

doesn't depend on the frequency, but only on the particle 
geometry, which is defined here by the parameters /^ and 
J. For imaginary part of cr( \\\, we have from Eq. (j48[) the 

linear enhancement with the frequency. This implies that 



-(1) 



30 I j^s 



J, 



(1)' 



(51) 



doesn't depend on the frequency, provided that the con- 
ditions uj <^ Vs and v <^ Vs are fulfilled. 



B. Size approach 

It remains to examine the cases of different relations 
between the particle sizes and the conduction electron 
mean- free path inside a particle. Using Eq. ((28|) . we can 
rewrite Eq. ((22|) in somewhat another form 



q^qi{e)^iq2i9), 



with 



QiiO) = 



{2k) ^o^'^+lw) «i 



(52) 



(53) 



sin2 6' 



where 2i?|| is the length of the MN along the z axis, which 
is directed along the principal spheroid axis, 2R± is the 
MN size along the x or y directions, and 



l^^ 



(54) 



has the sense of the length of an electron mean-free path 
(for Au at 0°C, e.g., I ~ 410A), and 



q2{0) = 



1 



#?) 



'sin^e 



' cos^ 6 



Here 



Vp Vp 

"" 2i?||' '^ 2Ri_ 



(55) 



(56) 



are the frequencies of the electron collision with the par- 
ticle surfaces along and across the 2;-axis of a spheroid, 
correspondingly. 

Let us consider the possible relations between / and 
particle sizes 2_R||, 2R±. 

i) The conduction electron mean-free path is much 
less than the sizes of the particle along particular direc- 
tions: I < 2i?j_, I < 2i?||. As follows from Eq. ((55t . 
qi{0)\v=vp ^ 1- In this case, an electron is scattered 
predominately inside the volume of MN. If, moreover, 
(72 — > 0, i.e., i^^n , Vs^± ^ to, then from Eq. ([?!]) . one gets 



* 



(57) 



Substituting Eq. ^7^ into Eqs. dM]), ^, and using 
Eq. (j331): we obtain the Drude-Sommerfeld formulas for 
real and imaginary parts of cr, presented above by the 
first term in the parenthesis of Eqs. (HI]), ((31]). 

ii) The mean-free path of a conduction electron is much 
greater than the particle size along particular directions: 
I ^ 2i?j_, I 3> 2i?||. In this case, an electron scattering 
occurs mainly from the inner surface of the MN. The 
electrons oscillate between the walls of the particle with 
different frequencies, excluding Vs,\\i ^s,x- In accordance 
with Eq. ([55]) . the inequality qi{0) <C 1 holds only for 
qi. The parameter 52 remains arbitrary. Formally, we 
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can put qi{9) — >■ 0, and for a real and imaginary parts of 
^-function one finds 

4 4/ 1 \ 

SR *((7)|gi^o = :j + 72 ( cosg2 sin 172 I , (58) 



3 g| 



and 



2 4 4/1 \ 

9 *((j)|gi^o = tH — 2 — cosg2 +sing2 • (59) 

12 92 92 \92 / 

The equations ([55]) , (|59l) allow to fulfil the calculation of 
the real and the imaginary parts of the ratio "ii / (v — iuj) , 
as we have done it before. Then, we obtain 
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-(1 - cos 92) 



(60) 
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<3l^0 



1 



^i 



cos 92 



1 



■ sm 92 



1 

3 9? \ ' 92 

(61) 

The parameter 92, defined by Eq. (|55l) . is governed by 
the frequency. Depending on the ratio between the inci- 
dent frequency and the frequencies h's,\\, J^s,±, the value 
of 92 can be greater or less than 1. This makes it difficult 
to perform subsequent analytical calculations of the in- 
tegrals involved into Eqs. ([M]) . ([55)) . Below, we dwell on 
some particular cases for which the calculations of a are 
the most simple. It should be noted also that the cor- 
responding expressions for the components of e can be 
easily obtained by substituting of cr into Eqs. (f32|) . ([33|). 



1. Conductivity of a spherical MN 

In the case of a spherical MN there are three char- 
acteristic frequencies which are considered usually: the 
frequency of an incident electromagnetic field w, the col- 
lision frequency of electrons in the particle volume v, and 
the vibration frequency between the particle walls Vs (if 
the particle size is less than the electron mean free path) . 
When ly > Vs, the mechanism of an electron scattering in 
the bulk is dominated, and an electron scattering from 
the particle surface gives only small corrections of the or- 
der of Vs/v. But we are interested in the case, when the 
mechanism of the surface electron scattering dominates, 
which corresponds to v < Vs- 

For particles of a spherical shape, the electric con- 
ductivity becomes a scalar quantity, and one can put 
i?ll — R±^ = i? in Eq. P7)) . then q and \l/-function in 
Eqs. (|25|) . ((26| are not depended on the angle 6. With 
an account for Eq. (|43| . this makes it possible to obtain 
for a the simple expression 



<^sph = ^11 = ^ I = 



Sne^ *(q) 
4 m V — iu) 



(62) 



with 9 = 2R{v — iu])/vF- To calculate o^^, one can use 



for different limit cases considered above. For example, 
we restrict ourselves here only to the case, when u <^ Vs- 
In this limit case, to a first approximation, one can put 
91 -^ 0. Then Eqs. ^^, ^Q with 92 = ui/vs can be 
used in Eq. ((62|) . As a result, one obtains for real and 
imaginary parts of a the following expressions 
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, (64) 



provided that v -^ Vs, where Vs = vp/{2R). The last 
expression in the Drude-Sommerfeld approximation looks 
like" 



sph 



^pi 



4tt v^ 



(65) 



If one puts the oscillation terms in Eqs. (|63|) . (|64|) equal 
to zero, and uses Eqs. ([5^ . (P5)) . then one obtains the 
expression for real part of the dielectric function which 
coincides with Eq. ([ij at i^ — >■ 0; but for imaginary part 
of the dielectric function, one gets Eq. ^ only if the 
replacement v — > ?iVs/2 will be done. This replacement 
is the same one as presented by Eq. ([3]), which has been 
often used in a phenomenological approximation. 




either Eqs. (gll), (gO]), or Eqs. (05]) , dH]) , or Eqs. ^, ^ 



oi 0.2 - 



FIG. 3. The real a' (solid lines) and imaginary cr" (dashed lines) 
parts of the ratio of an electric conductivity to the statical con- 
ductivity erg vs frequencies ratio, for a spherical Au particle with 
R = 200A. The thin lines correspond to the results obtained us- 
ing the kinetic method and the bold lines - the Drude-Sommerfeld 
formulas. 

The calculated results of the real and imaginary parts 
of an electric conductivity across the normalized to Vg 
frequency are shown in Fig. 3 for spherical Au nanopar- 
ticle, obtained with the use of both the kinetic method 
and the Drude-Sommerfeld formulas. The conductiv- 
ity is measured on the scale of the statical conductiv- 
ity a^ = ne^ /{mv). For calculations of cr, Eqs. (|34p . 



(|35| have been used in the kinetic case, and Eqs. ([2]) 
and (p5)) - in the Drude-Zommerfeld case. For illus- 
tration, the numerical parameters for Au particle^ and 
uipi = 1.37 X lO^^s"^ were chosen. As it can be seen 
in Fig. 3, the kinetic method appreciably changes the 
frequency dependence of asph at low frequencies, and 
at high frequencies (a; 3> Vs) gives the same result for 
asph as that which follows from the Drude-Zommerfeld 
formula. The difference between kinetic and Drude- 
Zommerfeld results is enhanced markedly as the parti- 
cle radius is decreased. The real part of a is peaked at 
uj/vs — >■ in both cases, whereas the imaginary part of 
CT - at w = I's in the Drude-Zommerfeld case, and at 
u! « 4i/s - using the kinetic method. 

For an extremely low to <^ Vs or an extremely high 
oj ^ Vs frequencies, one can find from Eqs. (p5|) . ([M)) 
after some algebra, the next simple approximation for a' 



sph 



and for ct'' 



sph 



167r 



^pi 



47r 






— ) , uj <t:vs 



(66) 



(67) 



The results of Eqs. ([55)) . (I57|) at uj :s> Vs correspond to 
those presented above by first terms in the square brack- 
ets of Eq. ([63t , (|64)) , accordingly. 

The calculations with employing of Eqs. (|63| . (|64l) for 
Au particle with R — 200A give the results similar qual- 
itatively to the ones above presented, but quantitatively 
they are at the a maximum of approximately 30% higher. 

In spite of the oscillation terms in square brackets of 
Eqs. dMl, dMl, the ratio of both cripjcrn and o'l^Ja^ 
does't oscillate with the frequency owing to the cut-off 
factors before these brackets. 
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FIG. 4. The dependence of e' for the spherical Au particle versus 
a radius R at the frequencies of lj =: 5.7 X lO^^s"^ (dashed line), 
6 X 10^''s-l (soUd Une), and 6.3 X Wi^^s'^ (doted Une). 

Now, let US discuss shortly the dependence of the di- 
electric permeability on the size of MN. In Fig. 4, we 



present the results of our numerical calculations of t(K) 
using Eqs. ([5^ . ([55]) at the fixed to. For illustration, 
we choose such frequencies from the frequency scale for 
which the above dependencies are the most pronounced 
for Au particle. 
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FIG. 5. The dependence of e" for the spherical Au particle across 
the radius R at the frequencies oi ui =: 5 X 10^*s^^ (dashed line), 
6 x 10l*s~i (solid line), and 7 X lO^^s-i (doted line). 

The real e' as well as the imaginary e" parts of the di- 
electric permeability oscillate, when the particle radius is 
increased. These oscillations have a damping character 
and practically are vanished for MNs of high radiuses. 
The real part of e tends to -1-1, and the imaginary part 
tends to at w — >■ 0. The real part has a first main 
minimum and the imaginary part has a first main max- 
imum at small values of R. Both the minimum and the 
maximum of e are slightly shifted towards the greater 
R with the frequency decreasing. In other words, this 
means that the resonance energy peak shifts toward the 
red-side with increasing sizes of Au particle. The oscilla- 
tion period in accordance with Eqs. (|55| . (IM)) is defined 
by 



r = 7r 



Vf 

R^ 



(68) 



and, as one can see, essentially depends on the product 
of Roj. It becomes shorter at Ruj <^ vp, and at Roj » vp 
- extends. Two types of oscillations may occur: at fixed 
Lu with varying of i?, or at fixed R with changing of u. 
For every fixed frequency there is the constant of e' or 
e" around which those quantities oscillate with altering 
of R. At high frequencies, the period of oscillations is 
sharply decreased and their amplitude are considerably 
lowered proportionally to the factors (ujpi/uj)'^ or ^^i/i^ 
before square brackets in Eqs. (j63p . (j64p . correspondingly. 
If we assume that the dielectric function of a surround- 
ing media e„i is close to unity, then for a spherical Au 
particle at a plasmon frequency w = ujpi/^/l + 2e,„ we get 
the closely packed oscillations of e'{R) within the ampli- 
tude interval —2 -;- 1. The behavior of <7'{R) and cr"{R) 
at the frequencies uj <^ Vg and lo ^ Vg can be seen from 
Eqs. dnH), dUl) as well. 
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2. Conductivity of an oblate MN 

For oblate particles 

/ > 2i?_L > 2i?||. 



(69) 



We shall consider, for convenience, the frequency inter- 
vals oj <^ Vs,i. and w ^ v^^^^. According to Eq. ([55]). 
the former interval corresponds to q2\v=VF < li and the 
latter one - to q2\v=VF > 1j respectively. It is easy to 
show that Eqs. (|60l) . (|6T|) for these frequency intervals 
transform to 



w < J^s I , 




(70) 



(71) 



Then, the integrals in Eqs. ([M|) . ([55]) can be calculated 
exactly in the approximations of ([70]). ([TT]) . and for an 
arbitrary spheroid aspect ratio between i?|| and R± we 
obtain the following expressions for the components of 
an electric conductivity 



/ / \ ^ pi 



U} <C l's,_L 
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(72) 
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(l)^ ^ 47r w 
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a; <C Vs.i^ 



1, w » I's 



(73) 



In the case of w <C i^s,±, Eq. ([7^ coincides with the first 
term found previously in Eq. (j47p . and in the case of 
u! ^ i/s ||, it coincides with the second term of Eq. (|4T|) 
at :^ = 0. Similarly, Eq. (1751) in the case of oj ^ z^^^ 
agrees with the first term found previously in Eq. (j48p 
with an accuracy of a numerical coefhcient, and in the 
case of oj :^ i^^ n coincides with the first term of Eq. ((42]) 
at z/ = 0. 

The above expressions can be transformed to their sim- 
plest forms for strongly deformed particles. Thus, for 
strongly oblate MNs, when R± ^ _R||, using Eqs. (jA3l) . 
(|A9|) from Appendix A, it can be easily found that at low 
frequencies 



9 , ,2 ^11 



pi Vf 



l^2f^-i 



w < iys,±, 



(74) 



and at high frequencies: 



647r \ uj ) J?,|| ' 



a'/2 



W > Z^5j|. 



(75) 



Similarly, for a" in the case of strongly oblate MNs at 
low frequencies, one gets 



2 '^11 V«ii 



^piR\\ 



arctan ■ 



i:j<z^s,_l, (76) 



and at high frequencies: 



cr_L = cr|| = 



47raj 



w » i/„ 



(77) 



It remains to consider, how the optical properties of 
MNs evolve at Vs,± < "^ < Vs,\\ between the low- and 
high-frequency interval. This interval is just that for 
which the parameter (72 can be greater or less than 1. 
In this case, only the numerical evaluations of the inte- 
grals entered into Eqs. (|M]). ((55)) can be performed. The 
obtained results for real part of a can be found in Ref. 
El. 



3. Conductivity of a prolate MN 



For prolate particles 



I > 2Rn > 2Ra 



(78) 



It is also worth noting that the approximations ((7D() . ((7T|) 
still remain to be valid for the case of inequalities (1751) , if 
the transposition Vg^^ ^ i^^ || is carried out in Eqs. ((70p . 
((7T|) . As a result, the integrals in Eqs. ((51]), ((5S]) can 
be easily calculated, and for prolate particles with an 
arbitrary aspect ratio of R\\/ R±, one obtains the results 
similar to ones given by Eqs. (l72l) . ((73]), in which the 
both replacements ((SO)) and ((A8() have been made for J 
and the integrals /, /^. Below, we write down only the 
results for a strongly prolate particles { R^ ^ R±) at low 
and high frequencies. Using Eqs. ((A4I) . (IA10[) . it is easy 
to find that 



rr' - -9-f.j2 R± 



aU2 



_ _9_ /"ifpiN^ vf_ 
~ 256 V w J B.± 

a'^ = 3cr[|/2 



UJ ^ Vs 



W > l's,±- 



(79) 



(80) 



For cr", just as for a' , in the case of strongly prolate 
MNs, one finds 



-L 167r y Vf I ' 



-2a'Ml + ln|^ 



a; ^ i/. 



(81) 



and at high frequencies a; ^ J^s,±, we have found a',! ~ 
(t" — uj'^i / {Attlu) , which coincides exactly with Eq. ((77)) 
for oblate particles. 
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Comparing Eqs. ([79| - ((8T|) for a strongly prolate par- 
ticle to the corresponding ones ([7^ - ([75)) for a strongly 
oblate particle shows that the character of the frequency 
dependence of the corresponding components of the elec- 
tric conductivity tensor remains practically the same. 
The behavior of this dependence reminds the asymptotic 
behavior of the Drude-Zommerfeld frequency dependence 
for an electron scattering in the volume of MN; this can 
be easily verified from the first summand in the right- 
hand sides of Eqs. (41), (42) by setting oj <^ v or oj :^ v. 
The difference consists only in the fact that for an elec- 
tron scattering in the volume, the high-frequency con- 
ductivity close to the frequency uj ^ v goes smoothly 
to the saturation. However, for the strongly asymmetric 
MNs there exists an entire transitional region between 
the minimum and the maximum transit rates, where the 
frequency dependence of a appreciably differs from the 
volume case. 

In the case of the frequency interval z^sj| < uj < i's,±, 
the numerical evaluations of the integrals involved into 
Eqs. (p4|). ([35l) can be performed and the obtained results 
for a have been presented in Ref. |45l |. 



V. RESULTS AND DISCUSSION 



is the radius of the MN, the greater is the correction to 
the frequency shift. 

Below, we will discuss the size dependence of optical 
properties of MNs in more details and will illustrate some 
results obtained in above sections. 
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FIG. 6. The frequency dependence of the reduced dielectric con- 
stant of Au particles with different R, A: 50 (thick lines) and 150 
(thin lines). Full lines are for a prolate particle with Rj_/R[\ = 0.1, 
and dashed lines are for an oblate one with Rj_/Ri\ = 10. 



One of the most important cases takes place for the fre- 
quencies Lu ^ Vs, with Vs — t^F/(2i?) 3> I'. In terms of q- 
language this corresponds to the situation when | gi | <^ 1 
and |g2| 3> 1. In this case, one may use the expressions 
([58l) . (|59|) at gi — > to obtain the results for an asymp- 
totically high value of g2- As a result, we get for a real 
and imaginary parts of an electric conductivity of a single 
spheroidal MN with an arbitrary aspect ratio of R±/R\\ 
the following expressions 



^(l)H-i(^)^ 



(1) 



{u) 



A-kKuj) 



9fF, 



9VF J 

4 R (!) 



(82) 



(83) 



provided that a; S> I's S> i'. The parameters /|| , I± are 
given in Appendix A by Eqs. (KJ) . (TMj) . 

For a MN of a spherical shape these expressions take 
the simplest form 









47r V w / 



(84) 



(85) 



provided that to ^ Vs ^ v. The last formula demon- 
strates simply that the electron interaction with the sur- 
face of a spherical MN can shift the imaginary part of a 
towards the red side of the frequency scale. The smaller 



In Fig. 6, the results of numerical calculations of the 
ratio between transversal and longitudinal components 
of the dielectric permeability versus the frequency are 
given for Au particles of different radiuses. The calcu- 
lations were performed with the use of Eqs. (|5^ . ([55)1 
and the same numerical parameters as above used. At 
high frequencies {uj ^ i^^), the ratio of ej_/eM — > 1 for 
different radiuses of Au particle. This result follows as 
well from Eqs. ([5^. ([75)1 a.t uj ^ Vg for particles of an 
oblate shape, and from Eqs. ((32)) . (|77)) - for particles of a 
prolate shape. At low frequencies (w <C Vg), the particle 
size is strongly effected on the ratio of i?e(e^)/i?e(e||). 
As it can be seen in Fig. 6, the curves for MNs of prolate 
and oblate shapes are merged into one with frequency 
grows, but for particles of a greater radiuses it becomes 
at smaller frequencies. 

The behavior of e' as a function of the ellipsoid aspect 
ratio R±/Rn at the frequency of a plasmon resonance in 
Au nanoparticle, embedded in the dielectric media with 
Em = 1, is plotted in Fig. 7. The numerical calculations 
with employing of Eqs. p2p . (I35p were performed for two 
radiuses of a nanoparticle and for two light polarizations. 
As one can see in this figure, the longitudinal component 
of e' in an oblate MN depends more stronger on the ratio 
R±/R\\ than the transverse one. It becomes especially 
pronounced for Au nanoparticle (with R = 50 A) of small 
oblateness (e.g., R±/R\\ < 15), when keeping the laser 
frequency fixed at uj — ujpi/y'l + 2em'- the oscillation 
magnitude of e^ (round the constant of e' = —2) becomes 
much larger than ones for ej,. The amplitude of these 
oscillations is enhanced and the period is extended with 
increasing of the ellipsoid aspect ratio R±/R^^. For MNs 
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of greater radiuses, the number of oscillations is decreased 
and their period is extended (compare the thin and thick 
curves in the Fig. 7 for Au particles with i? = 50 and 
150A, for example). The both components of e' no longer 
oscillate at high oblateness {R^/R\\ > 80) of Au particle 
with R ^ 50A. 
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FIG. 7. The dependence of the real part of the dielectric constant of 
Au particles with different i?, A: 50 (thick lines) and 150 (thin lines) 
on the ellipsoid aspect ratio, at a frequency of plasmon resonance 



oj = ajp;/V3~ 7.91 X 10^ 
for ±-component. 



^ . Solid line is for ||- and dashed line 



When the size of the particle is large enough, the ||- and 
_L-components of e' start to come together (see Fig. 2). 
For instance, at the frequency a; w 6 x lO^^s^^ it takes 
place for Au particle with R = 130A. At smaller frequen- 
cies it occurs for greater radiuses. 

Presented dependencies display mainly the behavior of 
e' for an oblate nanoparticle (i?_L/i?|| > 1). In the case 
of a prolate nanoparticle (0 < R^lR\^ < 1), the numeri- 
cal calculations with employing the same Eqs. (|32|) . (|35|) 
for Au nanoparticle of a small size (~ 50A) at the fre- 
quency of the plasmon resonance give the oscillations of 
a transversal component of e', the amplitude of which is 
enhanced and the period is reduced as soon as the pro- 
lateness of the MN is increased. The longitudinal com- 
ponent of eji oscillates in the prolate MN as well, but its 
amplitude is considerably smaller and the period is much 
larger than the proper ones for e'^. 

The imaginary part of the dielectric function as a func- 
tion of ellipsoid aspect ratio R±/R^^ is shown in Fig. 8. 
The numerical calculations were performed with the use 
of Eqs. p3l) . ([M)) at the plasmon frequency for two ra- 
diuses of nanoparticle and for two light polarizations. 
Since for prolate MN, the value of e" practically doesn't 
depend on the spheroid aspect ratio R±/R\\ (except for 
the case of a very high prolateness) , in Fig. 8 we present 
only the results for an oblate MN. In contrast to the size 
behavior of e', described above, one can see that the weak 
oscillations of e" hold together with linear increasing of 
e" just at small values of the aspect ratio R±/R\\ . For Au 
particles with R = 50A it takes place until R±/R\\ < 150, 



and is more sensible for the longitudinal component of 
e" than for the transversal one. The both components 
of e" reach the same maximum at some aspect ratios of 
Rj_/R\i, which value depends on the radius of particle. 
In the example depicting in Fig. 8 (for Au particle with 
R — 50 A), the values of ejf = e" ~ 1.8 have peaks at 
R±/R\\ > 140 for ||-component and at Rj_/R\\ ~ 700 
for _L-component of e". The another interesting feature 
of the dependence of e" on R^/ R\\ is the intersection of 
curves for ||- and _L-components for an oblate particle, as 
it takes place usually for the particle of a spherical shape. 




RJR« 

FIG. 8. The dependence of the imaginary part of the dielectric 
constant of Au particles with different R, A: 50 (thick lines) and 
150 (thin lines) on the ellipsoid aspect ratio, at the frequency of a 
plasmon resonance uj ~ 7.91 X 10^^ s^^ . Solid line is for ||- and 
dashed line for ^-component. 

Finally, we illustrate, how one may estimate the life- 
time of surface plasmon excitation (or any others) in a 
MN using the above derived formulas for a. As follows 
from our previous calculations2ii2S, the linewidth can be 
represented as 



r^(a;) = AiiLaaapiio), 



(86) 



where La is defined above after Eq. ([6]) and Uap is the 
real part of the conductivity tensor, given by Eq. (1231) for 
most general situations or by Eq. (p4)) for the spheroidal 
MN. For other particular cases, one may use for Re a 
Eqs. gl]), gll), dMl), dMl), dZa, dUl), (123, (I79|), (ini), and 
(|82p . presented above. In particular, considering only 
the nanoparticles of a spherical shape (L — 1/3) and 
restricting ourselves to the case i^ ^ :^s, we can choose 
Eq. ([55)1 for illustration. Substituting it into Eq. ([55]) . we 
obtain 



Iv, 




u 




W2 




sm 
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UJ 




Vs 





COS- 



^<")==S(?)' 

(87) 
Taking into account only the first term in (1871) , we recover 
the well- known ^— '•^^'•^^ 1/R dependence of F. 
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r.,.,«,^i^(^)^ 



As seen from Eqs. dST]), (|HH1), the lifetime (1/r) of an 
excitation in the MN depends not only on the nanopar- 
ticle radius, but also on the frequency (at which a given 
excitement is reasonable). For the frequency, which cor- 
responds to the excitation of surface plasmon in MN in 
a vacuum, uj = Wpj/yS, the following relation can be 
obtained from Eq. (j88p in energy units: 



Tl^{R) 



4 R 



(89) 



The oscillating terms in Eq. (I87t give rise to the oscil- 
lation of r around of Fq as a function of both the particle 
radius and the frequency. They can be represented at the 
frequency of surface plasmon as follows 



^SP 



(R) 



3a/3 h fvpY 

^^ V"rJ 

V3' 



. 2RuJr,l V'i Vp 

- sm ^ -I- 



V3 



VF 



2RUr, 



1 



2i?w„ 



V3 



Vp 



.(90) 



The amplitude and period of oscillations can be esti- 
mated by the values 



iVShvl 



-pmax _ 



T 



V3 



irvp 



Rcopi 



(91) 



accordingly. 

It is important to note that in kinetic method this os- 
cillatory behavior of F follows solely from the conditions 
of an electron scattering on the nanoparticle surface. 

Figure 9 shows the fuh linewidth F = T^^+T^f^, which 
is obtained by numerical evaluating of Eqs. ([M|) . ([TO)) for 
Na nanoparticles with parameter; 



iM. 



= 9.18 X 10 
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s~^ and vp = 1.07 x 10^ cm/s. One can see that oscillat- 
ing terms represent an important correction to Tq^{R) 
at small particle radii. This our result for Na nanopar- 
ticles only qualitative agrees with the similar results ob- 
tained in Refs. |37ll38| ). since we try to apply the kinetic 
method to the range of R, where the quantum effects 
(like, e.g., the Landau damping) play an important role. 
But, mostly, for the kinetics, the next inequality could 
be met 



i?> 



2TTh 

mvp 



(92) 



In order to study the significance of the oscillatory be- 
havior in more general situations, it is necessary to per- 
form the calculations of Eq. (1551) , using the real parts of 
Eq. (US]) or Eq. ([M]). This wiU be done separately. 




FIG. 9. Linewidth r(i?) of the surface plasmon resonance, as a 
function of radius for Na nanoparticles in units of Bohr radius ag ~ 
0.53A. The smooth term To{R) is given by Eq. I I89II (dashed) and 
the solid line corresponds to the sum of Eqs. I I89I1 and II90I I. 

There are several experimental data for a dielec- 
tric constant of powder of Ag and A\^ and for Ag 
nanoparticles42. For a single Au particle we have found 
only the experimental data for an optical response4^ A 
direct comparison of theoretical results with most of the 
available experimental measurements of the optical prop- 
erties of MNs are still a matter of debate because inho- 
mogeneous in nanoparticle size, shape and local environ- 
ment hide the homogeneous width of the surface plasmon 
resonance. There are very different data even for bulk 
permittivity of Au^i^ especially for its imaginary part. 



VI. CONCLUSIONS 

The kinetic equation method is used to study the pecu- 
liarities of the electron interactions with the surface of a 
spheroidal metal nanoparticle, when the electron scatter- 
ing from the particle surface becomes a dominant effect. 
The special attention was paid to study the modification 
of the Drude-Zommerfeld model applying to the optical 
properties of MN. The real and imaginary parts of the 
dielectric permeability at the frequencies above and be- 
low the characteristic frequency of a free electron passage 
between the walls of the particle were calculated for a sin- 
gle oblate and a prolate MN whose dimensions are much 
smaller than the wavelength of an electromagnetic wave. 

It was established for spherical MNs that the kinetic 
method appreciably changes the frequency dependence 
of electrical conductivity at low frequencies, and at high 
frequencies (w 3> Vg) gives the same result as one ob- 
tained from the Drude-Zommerfeld formula. Quantita- 
tively, the results obtained by kinetic method are of one 
or even two orders of magnitude lower (depending on the 
radius of MN) than those following from classical formu- 
las. The difference between results is enhanced markedly 
as the particle radius is decreased and the nanoparticle 
surface starts to play the more pronounced role. 



14 



The frequency dependencies of the components of the 
electric conductivity tensor a were found and their de- 
pendence on the spheroidal aspect ratio was investigated. 
Simple analytical expressions were found for this tensor 
in a strongly oblate or prolate MNs at low and high fre- 
quencies. 

The electron interaction with the surface of a spherical 
MN can shift the imaginary part of a towards the red side 
of the frequency scale. The smaller is the radius of the 
MN, the greater is the correction to the frequency shift. 

Two types of oscillations were established for small 
enough MNs: at fixed w with varying of i?, or at fixed 
R with changing of uj. These oscillations have a damp- 
ing character and practically are vanished else at high 
frequencies or for MNs of high radiuses. 

The ratio of /m(ej^)//m(e||) oscillates with increasing 
the frequency both for the prolate and the oblate MNs. 
In contrast, the ratio of the real parts of i?e(e_L)/i?e(e||) 
does not oscillate with frequency. Together with that, the 
real e as well as the imaginary e part of the dielectric per- 
meability oscillate, when the particle radius is increased. 
It was found that the particle size strongly effects the 
ratio of i?e(e^)/i?e(e||) at low frequencies (w ^ Vg)- 
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(A3) 



(A4) 



correspondingly, and in the limit case of R± = Rh — R, 
one gets 



(A5) 



We advance here as well the another typical integrals, 
which one meets under calculation of a. 
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Below, we present the values for integrals introduced 
in the Sec. IV. 



In the case of prolate particles {R\\ > R±), one should 
perform in Eqs. ([55]), (|il|) and in Eqs. (|XK|) . (|Mt the 
next replacement 
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For strongly oblate or prolate MNs {R± ^ R\\), using 
Eqs. (|A6|) . (jA7|) . it is easy to obtain that 
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- / ^sm^ev'{e)de 
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16i?|| 16i?|| i?i - Rf, 



X In 




(A2) 



For strongly oblate or prolate MNs, using the Eqs. 
, it is easy to find that 



(l),i?,i»i?,|i 2 R 1 ln2§^ - 1/2 



i),fl_L«-R|| ~ 4~R" ' \l/2 



(A9) 



(AlO) 



respectively, and in the limit case of R±^ = Rh = R, one 
finds from Eqs. (jM]) . (jM]) 



lim /,< = hm /^ = ^. 



(All) 
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